Average path length in uncorrelated random networks with hidden variables 
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Analytic solution for the average path length in a large class of uncorrelated random networks 
with hidden variables is found. We apply the approach to classical random graphs of Erdos and 
Renyi (ER), evolving networks introduced by Barabasi and Albert as well as random networks 
with asymptotic scale-free connectivity distributions characterized by an arbitrary scaling exponent 
a > 2. Our result for 2 < a < 3 shows that structural properties of asymptotic scale-free networks 
including numerous examples of real-world systems are even more intriguing then ultra-small world 
behavior noticed in pure scale-free structures and for large system sizes N — > oo there is a saturation 
effect for the average path length. 
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During the last few years random, evolving networks 
have become a very popular research domain among 
physicists 0, 0, 01 El- A lot of efforts were put into 
investigation of such systems, in order to recognize 
their structure and to analyze emerging complex prop- 
erties. It was observed that despite network diversity, 
most of real web-like systems share three prominent 
structural features: small average path length (APL), 
high clustering and scale-free (SF) degree distribution 
ULu>H,liLI9- Several network topology generators have 
beenproposed to embody the fundamental characteris- 
tics II Si IE El, Da. 

To find out how the small-world property (i.e. small 
APL) arises, the idea of shortcuts has been proposed 
by Watts and Strogatz [13]. To understand where the 
ubiquity of scale- free distributions in real networks comes 
from, the concept of evolving networks basing on prefer- 
ential attachment has been introduced by Barabasi and 
Albert ,6|. Recently Calderelli and coworkers ^3 have 
presented another mechanism that accounts for origins 
of power-law connectivity distributions. It is interest- 
ing that the mechanism is neither related to dynamical 
properties nor to preferential attachment. Caldarelli et 
al. have studied a simple static network model in which 
each vertex i has assigned a tag hi (fitness, hidden vari- 
able) randomly drawn from a fixed probability distribu- 
tion p(h). In their fitness model, edges are assigned to 
pairs of vertices with a given connection probability p^ , 
depending on the values of the tags hi and hj assigned 
at the edge end points. Similar models have been also 
analyzed in several other studies [T3. ITEl Il6| . 

A generalization of the above-mentioned network mod- 
els has been recently proposed by Boguna and Pastor- 
Satorras 0. In the cited paper, the authors have argued 
that such diverse networks like classical random graphs 
of Erdos and Renyi (ER), fitness model proposed by Cal- 
darelli et al. and even scale-free evolving networks intro- 
duced by Barabasi and Albert (BA) may be described by 
a common formalism. Boguna and Pastor-Satorras have 
derived analytical expressions for connectivity distribu- 
tions P(k) and relations describing degree correlations in 



such networks as functions of distributions of hidden vari- 
ables p(h) and the probability of an edge establishment 
Pij. In this paper we present an analytical description 
of main topological properties of the foregoing networks. 
We derive a general theoretical formalism describing met- 
ric features (i.e. APL, intervertex distance distribution) 
of random networks with hidden variables, assuming that 
the connection probability scales as p^ ~ hihj |l8j . The 
last assumption concerning the factorised form of p^ 
translates into the absence of two-point correlations and 
applies to a broad class of networks. 

The issue of the small-world property is of great im- 
portance for network studies. The property directly af- 
fects such crucial fields like information processing in dif- 
ferent communication systems (including the Internet) 
[T^. l2fl l2lL 123 . d isease or rumor transmission in social 
networks |23t l24l Eil as well as network designing and 
optimization |26t l27t I28L |29| . Not long ago, there was 
a strong belief that all the processes become more effi- 
cient when the mean distance between network sites is 
smaller. Recently however, it was shown that the small- 
world property may have an unfavorable influence on 
such phenomena like synchronizability |30| . 

Despite the universality and usefulness of the small- 
world concept, except a few cases [31], |33, |33l |34(, sat- 
isfactory calculations of the average path length (APL) 
almost do not exists. Even in the case of aged Erdos - 
Renyi graphs only a scaling relation (not an exact for- 
mula) describing APL is known Ier ~ In N/ ln(k) ||- In 
this paper we derive an exact formula for the average dis- 
tance Uj between any two nodes i and j characterized by 
given values of hidden variables hi and hj . Averaging the 
quantity Uj over all pairs of vertices we obtain the average 
path length characterizing the whole network. It is im- 
portant to stress that our formulas for APL do not posses 
any free parameters, therefore may be directly compared 
with results of computer simulations. In this paper we 
have tested our analytic results against numerical cal- 
culations performed for Erdos - Renyi classical random 
graphs, BA model and scale-free networks P(k) ~ k~ a 
with arbitrary scaling exponent a. In all the cases we 
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obtain a very good agreement between our theoretical 
predictions and results of numerical investigation. 
Let us start with the following lemma. 

Lemma 1 If A\ , A2, ■ ■ ■ , A n are mutually independent 
events and their probabilities fulfill relations VjP(-Ai) < e 
then 

n n 

P(|jA i ) = l-exp(-^P(A i ))-Q, (1) 



where 



o<0<E;ioM J /j!-(i + e)". 

Proof. Using the method of inclusion and exclusion [3£ 
we get 



p(U^) = £(-i) J ' +1 s(i), 

»=1 3=1 



(2) 



with 



S(j)= £ P{A il )P{A i2 )...P{A ij ) 

l<i 1 <i2<...<ij<n 



^(fl p ( A A -Qj. ( 3 ) 



where < Qj < (n?/j\ — (")) eK The term in bracket 
represents the total number of redundant components oc- 
curring in the last line of Neglecting Qj it is easy to 
see that (1— P(UAj)) corresponds to the first (n+1) terms 
in the MacLaurin expansion of exp(— P(j4j)). The ef- 
fect of higher-order terms in this expansion is smaller 
than R < (ne) n+1 /(n +1)!. It follows that the total er- 
ror of (JTJ may be estimated as Q < J2j=i Qj + This 
completes the proof. 

Let us notice that the terms Qj in J3J disappear when 
one approximates multiple sums X^i<» 1 <» 2 < 

by 

corresponding multiple integrals. For e = A/n -C 1 the 
error of the above assessment is less then A 2 exp(A)/n 
and may be dropped in the limit n — > 00. 

At the moment we briefly repeat (after Ref. ^tJ) the 
main properties of random networks with hidden vari- 
ables and connection probability pij given by 



Pij 



hihj 



(4) 



where (3 is a certain constant. In the case of random net- 
works, where two-point correlations at the level of hidden 
variables are absent we have 



(3 = (h)N, 



(5) 



whereas in correlated BA networks the prefactor gains 
another form. Boguha and Pastor-Satorras have shown 
that degree distribution P(k) in such uncorrelated net- 
works is given by 



p(*) = £ 



-P(h), 



(6) 



where p(h) describes a distribution of hidden variables. 
The above relation between both distributions P(k) and 
p(h) implies a relation between their moments 



(h n ) = (k(k-l)...(h-n + l)), 



and respectively 

(h) = (k), 



(h 2 ) = (k(k-l)). 



(7) 



(8) 



With respect to our following calculations the relation 
© requires a few comments. Firstly, let us note that 
for k — > 00 the Poisson-like propagator, that accompa- 
nies the distribution p(h) in the formula for P(k), is a 
sharply peaked function analogous to delta k- F° r this 
reason, in the limit of large nodes degrees we obtain a cor- 
respondence between the studied uncorrelated networks 
with hidden variables and random graphs with a given 
degree sequence (the so-called configuration model) j^il 



P(fc) ~ p(k). 



(9) 



Another very important conclusion that comes from 
considerations performed in Ref. 01 an d seems to affect 
our later derivations is that we can not generate uncorre- 
lated random networks with power-law degree distribu- 
tion P(k) ~ k~ a and the scaling exponent 2 < a < 3 by 
means of the factorised probability (@J (see also [sHIsS^. 
The axiomatic definition of probability requires pij < 1, 
giving the condition for the maximum value of the the 
hidden variable h max ~ v^V- When we think about hid- 
den variables as about desired degrees (as sketched in 
the previous paragraph) the condition for k max ~ h max 
is in contradiction to the cut-off of the power-law degree 
distribution k cut ~ N l /( a ~ l ) [3^ that allows for nodes 
with degrees higher than k max . For this reason, our for- 
malism describing metric properties of random uncorre- 
lated networks should not work well for SF networks with 
2 < a < 3. In contrast to the above discussion, we no- 
ticed that our analytical predictions are consistent with 
numerical calculations performed for scale-free networks 
with arbitrary scaling exponent a > 2. We suspect that 
the unexpected conformity for networks with 2 < a < 3 
may be related to the extreme small fraction of bad pairs 
of nodes with large degrees that do not fulfill the condi- 
tion p^ < 1 (see Appendix A). 

Now, we come back to the main subject of the paper, 
it means the issue of the average path length in random 
networks. Let us consider a walk of length x crossing 
index-linked vertices {i, Vi,V2 ■ ■ ■ vr x _i),j}. Because of 
the lack of correlations the probability of such a walk is 
described by the product Pi Vl Pviv 2 Pv 2 va ■ ■ -Pv^^j, where 
Pij gives a connection probability between vertices i and j 
Q . At this stage it is important to stress that the graph 
theory distinguishes a walk from a path [40|. A walk is 
just a sequence of vertices. The only condition for such a 
sequence is that two successive nodes must be the nearest 
neighbors. A walk is termed a path if all of its vertices 
are distinct. In fact, we are interested in the shortest 
paths. In order to do it, let us consider the situation when 
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there exists at least one walk of the length x between 
the vertices i and j. If the walk(s) is(are) the shortest 
path(s) i and j are exactly x-th neighbors otherwise they 
are closer neighbors. In terms of statistical ensemble of 
random graphs [ilj the probability pij (x) of at least one 
walk of the length x between i and j expresses also the 
probability that these nodes are neighbors of order not 
higher than x. Thus, the probability that i and j are 
exactly x-th neighbors is given by the difference 



Pijfr) =Pij{x) ~Pij(x - 1). 



(10) 



In order to write the formula for p^ (x) we take advan- 
tage of the lemma 



Pij{x) = 1 - exp[ 



N 
v\ — 1 



N 



2^ PiV! ■■■Pv (x ^ 1) j], (11) 



where N is the total number of vertices in a network. 
A sequence of (x + 1) vertices {i, v\, V2 ■ ■ ■ , be- 
ginning with i and ending with j corresponds to a sin- 
gle event Ai and the number of such events is given by 
n = N x ^ 1 . Putting @ into 111(1 and replacing the sum 
over nodes indexes by the sum over the hidden variable 
distribution p(h) one gets 



pij(x) = 1 - exp 



hihj f(h 2 )N 



(h 2 )N V (3 



(12) 



Due to l|10l) the probability that both vertices are exactly 
the x-th neighbors may be written as 



p* j (x) = F(x-l)-F(x), 



where 



F(x) 



exp 



(h 2 )N 



hi hi 



(h 2 )N V 



(13) 



(14) 



The above calculations require a few comments. First 
of all, note that the assumption underlying JQ) is the mu- 
tual independence of all contributing events Ai. In fact, 
since the same edge may participate in several x— walks 
there exist correlations between these events. Neverthe- 
less, it is easy to see that the fraction of correlated walks 
is negligible for short walks (x <C N) that play the major 
role in random graphs showing small- world behavior. It 
is also important to stress that our formalism does not 
neglect loops. 

Let us point out that having relations (1K'>[1 and (|14fl . 
describing the probability that the shortest distance be- 
tween two arbitrary nodes i and j equals x, one can 
find almost all metric properties of studied networks 
[42| . For example, averaging (|13fl over all pairs of ver- 
tices one obtains the intervertex distance distribution 
p(x) = ((p*j(x))i)j. It is also possible to calculate z x 
- the mean number of vertices a certain distance x away 
from a randomly chosen vertex i. The quantity can be 
written as z x = J p*j(x)p(hj)Ndhj. Note that taking 



only the first two terms of power series expansion of 
both exponential functions in (|13fl and making use of Q 
and JHJ one gets the relationship z x = z\{zij zxf^ 1 — 
(k)((k 2 ) / (k) — l) x_1 that was obtained by Newman et 
al. |36j when assuming a tree-like structure of random 
graphs with arbitrary degree distribution. 

Taking advantage of l|13|l one can calculate the expec- 
tation value for the average distance between i and j 



kjihuhj) = V xp*Ax) 



F[x). 



(15) 



Notice that a walk may cross the same node several times 
thus the largest possible walk length can be a; = oo. 
The Poisson summation formula allows us to simplify the 
above sum (see Appendix B) 



lij {hi , hj) 



\nhih 3 ; +ln7V + ln(/i 2 ) -7 1 
In N + ln(/i 2 ) - hi /3 + 2' 



(16) 



where 7 ~ 0.5772 is the Euler's constant. The average 
intervertex distance for the whole network depends on a 
specified distribution of hidden variables pih) 



I 



-2{hih) +lniV + ln(/i 2 ) _ 7 1 
\nN + \n(h 2 ) -ln/3 + 2' 



(17) 



We need to stress that both parameters Uj and I diverge 
when the argument of the logarithmic function in the de- 
nominator of both expressions (|16fl and (|17(l approaches 
one i.e. N(h 2 )/(3 = 1. Note, that substituting © for 
(3 in the last condition and then taking advantage of Q 
one recovers the well-known estimation for percolation 
threshold (k 2 ) / (k) = 2 in undirected random networks 
with arbitrary degree distribution pol l43l lil . l45j (see 
Appendix C). 

To test the formula (|17|) we start with the well-known 
networks: ER classical random graphs, BA model and 
scale-free networks. The choice of these networks is not 
accidental. The models play an important role in the 
network science. The ER model was historically the first 
one but it has been realized it is too random to describe 
real networks. The most striking discrepancy between 
ER model and real networks appears when comparing 
degree distributions. As mention at the beginning of the 
paper degree distributions follow a power-law in most 
of real systems, whereas classical random graphs exhibit 
Poisson degree distribution. It was found that the most 
generic mechanism driving real networks into scale-free 
structures is the linear preferential attachment. The sim- 
plest model that incorporates the rule of preferential at- 
tachment was introduced by Barabasi and Albert 0. 
Other interesting mechanisms leading to scale-free net- 
works were proposed by Goh et al. and Caldarelli et 
al. 01 • Goh and coworkers were the first who pointed 
out that power-law connectivity distribution P(/c) may 
result from the Zipf law applied to hidden variable dis- 
tribution pih) ~ h~ a . The concept of the Zipf law has 
been next developed by Caldarelli et al. in their pa- 
per |l2^. In fact, the most important achievement of the 
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FIG. 1: The average path length Ier versus network size N 
in ER classical random graphs with (k) = 4, 10, 20. The solid 
curves represent numerical prediction of Eg. 11911 . 



mention paper by Caldarelli et al. relates to a nontrivial 
discovery that scale-free networks may be also obtained 
from exponential distribution of fitnesses p(h) <~ e . 
Since however, the case of scale- free networks with ex- 
ponentially distributed fitnesses does not fulfill (@J, we 
do not take it into account in this paper. In the present 
study, we examine the case of scale-free networks with 
underlying scale-free distributions of hidden variables. 

Below we show that our formalism describing metric 
properties of random networks may be successfully ap- 
plied to all the above listed network models. 

Classical ER random graphs. Note that the only way 
to recover the Poisson degree distribution form the ex- 
pression 10 is to assume 



PEB,{h) = 5( k ) <h . 



(18) 



Now, applying the distribution PER.(h) to i|17|) we get the 
formula for the average path length in classical random 
graphs 



ER 



In N - 7 1 
ln(fe) + 2' 



(19) 



Until now only a rough estimation of the quantity has 
been known. One has expected that the mean intervertex 
distance of the whole ER graph scales with the number 
of nodes in the same way as the network diameter. We 
remind that the diameter d of a graph is defined as the 
maximal shortest distance between any pair of vertices 
and dER = lniV/ln(fc) [3j. Fig^ shows the prediction 
of the equation (|19|) in comparison to the numerically 
calculated APL in classical random graphs. 

Scale-free BA networks. The basis of the BA model 
is its construction procedure 0, |46| . Two important in- 
gredients of the procedure are: the continuous network 
growth and the preferential attachment. The network 
starts to grow from an initial cluster of m fully connected 
vertices. Each new node that is added to the network cre- 
ates m links that connect it to previously added nodes. 
The preferential attachment means that the probability 



of a new link growing out of a vertex i and ending up in 
a vertex j is given by 



(20) 



where kj(ti) denotes the connectivity of a node j at the 
time ti, when a new node i is added to the network. Tak- 
ing into account the time evolution of nodes degree in BA 
network (i.e. putting kj(ti) = m^J ti/tj), the probability 
of a link between i and j can be rewritten in the following 
form 



in 



1 



(21) 



that is equivalent to 
hj = 1/yfEj and j3 B A 



(0J when assuming hi = l/y^, 
= 2/m. The distribution of hidden 



variables pba (h) in BA networks follows the relation 



PBA{hi)dhi = P{ti)dti 



(22) 



where P(U) — 1/N is the distribution of nodes attach- 
ment times ti for a network of size N . After a simple 
algebra one gets 



PBA(h) = —h~ 



(23) 



for h — 1/yN, . . . , 1. Now, it is simple to calculate the 
average distance l(TT)|) between any two nodes in BA net- 
works 



hi(hihj) — ln(m/2) — 7 3 
lnlniV + ln(m/2) + 2' 



(24) 



Averaging l|24|l over all pairs of vertices one obtains APL 
characterizing the whole network 



BA 



IniV- ln(m/2) - 1 -7 3 
lnlnA + ln(m/2) + 2' 



(25) 



Fig|21 shows the APL in BA networks as a function of 
the network size N compared with the analytical formula 
B25JI . There is a visible discrepancy between the theory 
and numerical results when (k) — 2m = 4. The discrep- 
ancy disappears when the network becomes denser i.e. 
when (k) increases. The same effect will appear at Fig|U 
letting us deduce that for some reasons our formalism 
works better for denser networks. 

Scale-free networks with arbitrary scaling exponent. 
Let us start with the well-known model of scale-free net- 
works introduced by Goh et al. (Model A) and its 
certain modification proposed by Caldarelli et al. (Model 
B) [l2j. We show that both models A and B possess 
peculiar properties that make application of our theoret- 
ical approach impossible. Next, we make use of a general 
procedure described at the beginning of the paper to gen- 
erate uncorrelated networks with asymptotic power-law 
connectivity distributions (Model C). 
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N - network size k 



FIG. 2: Characteristic path length Iba versus network size N 
in BA networks. Solid lines represent Eg. 1251 . 



Model A. To construct the network one has to per- 
form the following steps: (i.) prepare a fixed number iV 
of vertices; (ii.) assign fitness (hidden variable) hi = i~ r , 
with < t < 1, to every node i = 1, . . . , N; (iii.) select 
two vertices i and j with probabilities equal to normalized 
hidden variables, hi/ ((h) N) and hj/((h)N), respectively, 
and add an edge between them unless one already exists; 
(iv.) repeat previous steps until mN edges are made in 
the system. Goh and coworkers have showed that the re- 
sulting network generated in accordance with the above 
procedure exhibits asymptotic power-law degree distri- 
bution 



where 



P(k) 



a = 1 + 1/7 



(26) 



(27) 



that gives 2 < a < 00. Although in these networks 
probability of a connection approximately factorizes (@J 



1 



1 



hihj 



■iN 



hihj 



(28) 



where (3 a = (h) 2 N/m, there is one important feature 
of the model. The non- analytic statement, included in 
the step (iii.) of the construction procedure expressed 
as add an edge unless one already exists, gives rise to 
uncontrolled intervertex correlations both for large m and 
small a < 3. 

Model B. Caldarelli and coworkers have modified the 
original model introduced by Goh et al. by assigning to 
nodes random fitnesses hi taken from a given distribution 
p(h) ~ h~ a , instead of deterministic values hi = i~ T . 
They also assumed a modified edge establishment pro- 
cess: for every pair of vertices i and j a link was drawn 
with probability where (3b = (hmax) 2 ■ Although the 
foregoing value of (3b assures us of < 1, it is strongly 
overestimated and makes resulting networks very sparse 
with a large content of isolated nodes |47|. 

Model C . In order to avoid features incorporated in 
both models A and B, we have generated networks pos- 
sessing asymptotic scale-free behavior for fc 3> 1 coming 



FIG. 3: Model C. Degree distribution P(k) (the main layer) 
for random network with underlying hidden variable distribu- 
tion given by a power-law (the inset). Scatter data represent 
results of numerical calculations, whereas solid curves express 
formulas © and l|29|l . respectively for P(k) and p(h). 



out of power-law distributions of hidden variables 



p(h) 



(a - l)m( a -V 



h a 



(29) 



for h = m, . . . , h max , where h max ^ mN 1 / 1 * - ^ (see |^) 
and connection probability given by (£Q) and JSJ. A typ- 
ical behavior of connectivity distribution P(k) for net- 
works generated in accordance with this procedure is pre- 
sented at Fig|3J Note that for fc > to the connectivity dis- 
tribution is well described by the power law P(k) <~ k~ a 

©■ 

Applying the distribution i(29"|) to the formula l(T7j) one 
obtains 



for a > 3 



;a>3 



lniV + hi(f- 



for a = 3 



ia— 3 



ln(!=§) + lnm 



lnJV-ln(f)-l-7 
lnlniV + ln(f ) 



1 

2' 



3 

2' 



(30) 



(31) 



• for 2 < a < 3 

2 



ja<3 = 



(^ T )lnJV + ln(g^)-(^ T )-7 , 1 
(|=f)lnJV + lii(fE5)+lnm 



2- ^) 



Fig0| shows predictions of the above equations in com- 
parison with numerically calculated shortest paths. We 
would like to stress that regardless of the value of a, for 
denser networks (with higher values of parameter to) , one 
can observe an excellent agreement between our theory 
and numerical results. 

Summarizing, depending on the value of scaling expo- 
nent a one can distinguish three scaling regions for the 
average path length in scale-free networks. In the limit 
of large systems N — ► 00 , APL scales with network size 
according to relations 
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FIG. 4: Model C. The average path length versus network 
size N for a = 4 (CI), a = 3 (C2) and a = 2.5 (C3). The 
scatter data represent numerical calculations. Solid curves 
with open squares in the case of m — 5 (open circles in the 
case of m = 20) express analytical predictions of Eqs. llMIH . 
JSP and respectively for (C1),(C2) and (C3). 



• for a > 3 



• for a = 3 



• for 2 < a < 3 



Z Q>3 - lniV, 



/a— 3 



;a<3 



lniV 
In In AT' 



2 1 



3-a 2 



(33) 



(34) 



(35) 



Note that although the results for a > 3 are consistent 
with estimations obtained by other authors [33l l3^ | , the 
case of 2 < a < 3 is different. In opposite to previous 
estimations signaling the double logarithmic dependence 
l a<3 ~ In In N, our calculations for the same range of a 
predict that there is a saturation effect for the mean path 
length in large networks. Since the assumption underly- 
ing estimations leading to double logarithmic dependence 
in APL was a pure scale-free behavior of degree distri- 
bution, we suspect that this discrepancy may result from 
ambiguous behavior of P(k) in our model. Let us note 
that in our model C there is a relatively small number of 
nodes with small degrees k (see Fig|3J). Since distances 
between such nodes are usually very large in compari- 
son to distanced between nodes with higher degrees, thus 
their absence may lead to the domination of the APL 
parameter by distances between the population of nodes 
characterized by medium degrees. Our result shows that 
for 2 < a < 3 structural properties of asymptotic scale- 
free networks including numerous examples of real- world 
networks may be even more intriguing then ultra-small 
world behavior reported for pure scale- free systems. 

To conclude, in this paper we have presented theoret- 
ical approach for metric properties of uncorrelated ran- 
dom networks with hidden variables. We have derived a 
formula for probability p*j (x) (|13|l that the shortest dis- 
tance between two arbitrary nodes i and j equals x. We 
have shown that given p*j(x) one can find every struc- 
tural characteristic of the studied networks. In partic- 
ular, we have applied our approach to calculate exact 
expression for the average path length (|17|l in such net- 
works. We have shown that our formalism may be suc- 
cessfully applied to such diverse networks like classical 
random graphs of Erdos and Renyi, evolving networks in- 
troduced by Barabasi and Albert as well as random net- 
works with asymptotic scale-free connectivity distribu- 
tions characterized by arbitrary scaling exponent a > 2. 
In all the studied cases we noticed a very good agree- 
ment between our theoretical predictions and results of 
numerical investigation. 

Acknowledgments. First of all, we wish to thank to 
anonymous Referee for critical comments which helped 
us to refine our paper. In the preliminary version of the 
paper [4Sj we have used a node degree notation, that was 
not accurate enough with reference to the analyzed prob- 
lem. Following the comments given by the anonymous 
Referee we have reformulated our approach in a language 
of hidden variables. We are also thankful to Sergei Doro- 
govtsev for his critical comments to an earlier draft of 
this paper. One of us (AF) acknowledge The State Com- 
mittee for Scientific Research in Poland for support under 
grant No. 2P03B01323. 

Appendix A. The condition pij < 1 @| is not fulfilled 
for pairs of vertices i and j possessing large hidden vari- 
ables (or desired degrees) hi and hj. To justify our cal- 
culations, we have to assure ourselves that the fraction 



7 



of such pairs is very small 



p(hj) / p(hi)dhidhj <C 1. 

min Jpij/3/hj 



(36) 



one realizes that in most of cases 



H ,0j 







Using the Chebyshev's inequality [35| and solving 136(1 
with respect to pij < 1 one gets 



(h 2 



(h 2 )N 

that gives F(0) = 1. One can also find that 



(40) 



(hy 



■((h 2 ) - (h) 2 ) « iV 2 , 



(37) 



where we assumed (3 — (h)N . ft can be shown that 
every network that is considered in this paper fulfill the 
condition. 

Appendix B. The Poisson summation formula states 



J F(x)dx = -Ei f- 



hlhA/ nr^r (4i) 



(h 2 )N 







oo n 

J2F(x) = -F(0) + 

poo °° / P oo 

J F(x)dx + 2 ^2 F ( x ) cos(2mrx)dx 



(38) 



Applying the formula to 115(1 

oo 

kj(hihj) = y^exp 



x=0 



/i^ ft -J 


f{h 2 )N\ x l 


(h 2 )N 


V (3 J \ 



(39) 



where Ei(y) is the exponential integral function that for 
negative arguments is given by Ei(—y) = 7 + bay [Ifij . 
where 7 ~ 0.5772 is the Euler's constant. Finally, it 
is easy to see that owing to the generalized mean value 
theorem every integral in the last term of the summation 
formula 1(38(1 is equal to zero. It follows that the equation 
for the APL between i and j is given by 1(16(1 . 

Appendix C. Note that, using additional assumptions 
one can simply reformulate both formulas l(16|) and l(17|) 
as well as the whole formalism in terms of node's degrees 
instead of hidden variables. For more details sec 481. 
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